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I. INTRODUCTION 

Cosmological observations put strong constraints on processes which could have occurred in 
the early universe. For example, models of inflation are tested with the properties of the cosmic 
microwave background radiation (CMB), such as the CMB anisotropies, CMB polarisation, non- 
Gaussianity and spectral distortions to the black-body spectrum. The Planck satellite mm 
provides the most recent observational data of the CMB. Non-Gaussian statistics originating from 
inflationary physics can furthermore be probed with studies of the large scale structures (LSS) in 
the universe |3]. This is one of the goals of state-of-the-art current and future experiments, such 
as the Dark Energy Survey [Ij, the Large Synoptic Survey Telescope (LSST) [5] and the Euclid 
satellite [B]. 

The number of different inflationary models is vast. The simplest models consist of a single scalar 
field minimally coupled to gravity. However, the phenomenology of even these simple models is 
rich, with hundreds of different choices for the inflationary potential [7]. More complex models 
arise from including more than one scalar field, which can lead to qualitative differences. These 
differences occur due to fluctuations not just in one field direction, but now in more than one 
direction (i.e. isocurvature or non-adiabatic pressure perturbations, to which we will return later). 
Then, one could consider a single field with a non-standard kinetic term, such as k-inflation [8]; 
these such models are often motivated by theories with extra dimensions, e.g., DBI inflation from 
string theory |9]. For more complexity, the scalar field could have a non-minimal coupling to 
gravity, such as in the newest version of the Higgs inflation model [lOj (however, for the majority 
of cases, this can be treated as a field with a minimal coupling by moving to the Einstein frame and 
modifying the potential). Finally, the most complex inflationary models contain multiple scalar 
fields non-minimally coupled to gravity, and with non-standard kinetic terms mm- 

By comparing each model’s predictions with observational data, we can rule out regions of 
model space, with the ultimate goal to obtain a single inflationary model which best fits the data. 
Recent data provides bonnds on the gravitational wave signature for which the simplest single 
field inflationary model with an m?4>^ potential is disfavoured m- Therefore, it is particularly 
important to continue to investigate the dynamics and observational predictions of inflationary 
models beyond the simplest single scalar field model. One interesting model not belonging to the 
single field class is the curvaton model mm- This model consists of a second field, the curvaton, 
in addition to the inflaton. The curvaton is dynamically unimportant during inflation, but its 
fluctuations source the curvature perturbation. 
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In this paper, we address the question of whether a possible stabilisation of the Planck mass 
(or gravitational constant) just after inflation can have a sizeable effect on the primordial power 
spectrum of the curvature perturbation. In theories in which the four-dimensional Planck mass 
are not constant, its dynamics is usually driven by one or several moduli fields. These describe 
for example the size of the extra-dimensional space. Since the time evolution of the gravitational 
constant is strongly constrained by experiments (see e.g. HZ! for a recent update on experimental 
tests of General Relativity), any stabilisation of the moduli field(s) must have happened in the early 
univers^ The stabilisation could have happened well before inflation ended, affecting scales well 
outside the visible horizon. If the stabilisation happened during the last 60 e-folds of inflation, 
possible signatures in the curvature perturbation power-spectra can be produced [T8|. If the 
stabilisation happened later, in the radiation dominated epoch, the rapid oscillations of the scalar 
field(s) can affect the expansion history |19l 120] . In the scenario discussed in this paper, Newton’s 
constant stabilised a few e-folds after inflation. We take into account the possibility that the 
moduli field driving the evolution of Newton’s constant can decay into other particles as well. Our 
setup is therefore related to the standard curvaton scenario. In the absence of a non-adiabatic 
pressure perturbation, the curvature perturbation ^ on uniform density hypersurfaces is known 
to be conserved on superhorizon scales [21]. However, this is not necessarily the case if several 
fields are dynamically important. Even if inflation has ended, the decay of fields at a later stage 
can significantly influence the evolution of C (see e.g. jHHISlISlHZS]). In the case of a scalar field 
driving the evolution of the Planck mass, we find that the non-minimal coupling to the Ricci scalar 
can boost the amplitude of the curvature perturbation by several orders of magnitude, even if the 
Planck mass varies only by a very small amount. 

The paper is structured as follows: in the next section, we present the model, and the governing 


equations for the background and perturbations of the model. Then, in Section III we describe our 
numerical procedure, before presenting results in Section IV Finally, we conclude in Section jV] 


II. THEORY AND ANALYTICAL RESULTS 


The model we consider consists of two scalar fields, namely of the inflaton (p and the field a, 
which describes the evolution of the Planck mass. In the Jordan frame, the action is given by 


S= 


^f{cr)R - + d^ad^a) - V (p, a) + Ant 


( 1 ) 


^ Alternatively, the post-inflation evolution of the field(s) could be very slow. 
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where is the metric tensor, V{(p, a) the potential and £int is the interaction Lagrangian, de¬ 
scribing the perturbative decay |26) of both the (p and a into radiation. By working in the Jordan 
frame, the decay rates can be calculated in the standard way. We denote them by and F”" 
respectively. Since we are interested in the effect of stabilising the Planck mass, we expand f{cr) 
around its minimum, keeping only the leading term: 


a 


f{(^) = 1 + 


( 2 ) 


We denote the masses of the fields by and and assume that the fields are not directly 
interacting. Note that we are working in units with reduced Planck mass Mpl = 1.Therefore, the 
potential is given by 


^O') = \rnl(p‘^ + 


(3) 


To consider the evolution of cosmological perturbations, we work in the longitudinal gauge, in 
which the line element takes the form |27ff29] 

(4) 


ds^ = —(1 -b 2^)dt^ + a^{t) (1 — 2T) dijdx^dx^ 


Here, a{t) is the scale factor, $ and T are independent metric perturbations, which depend on all 
coordinates. In the Jordan frame, $ and T are not equal even in the absence of anisotropic stress 


(see Eq. (15) below). The equations of motion for the system can be obtained by varying the 
action in Eq. (0. In the background, we have evolution equations for the two scalar fields 

ip = -V^-2,H^-T^^, (5) 

a = -Va- 3H& + RU/2 - F'"d , (6) 


in addition to an energy conservation equation for the radiation fluid 

= -AHp^ + + F'^d^ , 


and the Friedmann equation 


^ 3/ 


a2 ^2 


y + y TE + Pt, 


fgdH 

f 


(7) 


( 8 ) 


In these equations, a subscript cp or a denotes a partial derivative with respect to the field, and 
we have written the derivatives with respect to cosmic time, t. We shall also use the slow roll 
parameter defined by m, 

H 




H2 ’ 


(9) 
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in order to more simply write the Ricci scalar, which can be expressed as R = 6H^{2 — e) along 
with its perturbation 


h‘2 

SR = + 4^') - 2R4> + 2^(4- - 24') 


( 10 ) 


Considering now the linear perturbations, we obtain a Klein-Gordon equation for each field m 


(11) 


Scj) = - ^ + R00 + <(.(4> + 34-) - 


a2 


faaR 


5a = -3H5a - ( ^ + ) <5a - V^^5^ + + 3’f) - 2R,$ + (2R4> + 5R ), (12) 


along with a conservation equation for the radiation fluid 




6p^ = -AH5p^ + 4p^4> - 2^(^ + HA>) + 2r^(,/.5(/- - ^$) + 2r^{a5a - ^$). (13) 

The metric perturbation 4^ satisfies the following evolution equation 
T = -3F4' H^{3 - 2e)4> 

<j)5(p + a6a - + a^)A> - V^Scj) - VJa - 2/$ - /(T + 2F4>) 


1 


-J (^ + ^-V + f + 2Hf] +Sf + 2H5'f+’^5f 


(14) 


along with the constraint 


4> = 4' - 


SI 

f 


(15) 


The predictions from inflationary models are mapped onto observations (such as the temperature 
fluctuations of the CMB) in a simple way by introducing a curvature perturbation. The curvature 
perturbation on uniform density hypersurfaces, is defined as 


C = -^--Sp, 
p 


(16) 


where here, p and 5p are the energy density and perturbation for the entire matter content of the 
universe. We can obtain an evolution equation for C, which, in the large-scale limit, takes the form 

H 


c = -- 


:5Pna.d , 


(P + P) 

where the non-adiabatic pressure perturbation, (5Pnad) is defined as 

^Pnad ^5P - Ksp. 

P 


(17) 


( 18 ) 
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Table I: A table clarifying our notation for the subscripts denoting various stages in the evolution 

of a. 


dini 

The initial value of a 

<^end 

The value cr reaches at the end of inflation, before rolling down to its minimum 

^min 

The minimum in the expansion of /(cr) 


For a minimally coupled single field model of inflation (or, equivalently, for a universe containing 
a single fluid), the curvature perturbation, is conserved for both canonical and non-canonical 
models of inflation, independent of the theory of gravity [ZU Ea Eg. This allows us to compare 
inflationary predictions directly to observations by mapping the field fluctuations onto the curvature 
perturbation. Since it is conserved, we do not need to worry about the mechanism by which inflation 
ends and the universe reheats. 

However, moving beyond these simple models, the non-adiabatic pressure (or entropy) perturba¬ 
tion is non-zero, and therefore the curvature perturbation can continue to evolve and be enhanced 
on super-horizon scales. This feature has been exploited in numerous scenarios containing multiple 
minimally coupled scalar fields (see, e.g.. Refs. [Z3HM1IMMSS] ) • In these models, we must take into 
account the reheating phase in order to make reliable predictions. Models with non-minimally cou¬ 
pled scalar fields, on the other hand, will produce a distinct source of entropy perturbations, arising 
due to the coupling of the scalars. In the model we present above, it is expected that, after inflation 
has ended and during the reheating phase, these entropy perturbations can become sizeable due 
to the fact that / and / no longer need to remain small m- It is these non-adiabatic pressure 
perturbations, and the resulting amplification of the power spectrum, that we will investigate in 
the remainder of the paper. 

We will solve the system of equations derived above. It will be necessary to follow the evolution 
of the secondary field, a, throughout the inflationary phase, the decay of the inflationary field, 
(j), and the radiation epoch right up until a itself has decayed and no longer contributes to the 
overall energy density of the universe. This is important so as not to restrict ourselves to the case 
where the auxiliary field starts at its minimum, cimin (see Table [I] for subscript notation), and is 
pushed away, but to also include cases where the field evolves towards cimin during inflation. It 
is often the case that when studying subdominant, curvaton-like fields, the calculation begins in 
a post-inflationary radiation-dominated phase and proceeds from there. For the usual, minimally 
coupled fields this is sufficient, since a does not evolve until late on, after the end of inflation |15] . 
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However, as this no longer holds in our case, we must track it throughout. It is still important 
for this coupling to remain small so as not allow a to contribute too much and impact upon the 
dynamics of inflation itself. 

This feature, which distinguishes our setup from the standard curvaton scenario, arises due to 
the explicit coupling to the Ricci scalar causes the field to obtain an effective mass, which might 
not be small compared to H^. In the slow-roll approximation we find that cr evolves according to 


3ai=^N 
a (X e 3-'= , 


(19) 


where e is assumed to be roughly constant and N is the e-fold number. This equation follows 
directly from Eq. Q, writing this equation in terms of e-fold number N and neglecting the bare 
mass of the field. Therefore, it is often the case that the field-value of a at horizon crossing is 
different from the value of a at the end of inflation. 


III. NUMERICAL METHOD 

We will solve the governing equations derived in the previous section using a code written in 
Python. This starts at the beginning of inflation and runs right through to the end of the decay 
of the second field, when the power spectrum has reached its final value. This is the full numerical 
simulation including field perturbations, their gravitational counterparts and those in the radiation 
fluid created in the final stages. We begin by running through the background in order to ascertain 
the values needed to set up the initial conditions for each mode. A:, such that /c* = 50aiHi. We then 
set the initial conditions of the perturbations as those of the Bunch-Davies vacuum m at this 
point and begin the full perturbation code. The perturbation equations are then each integrated 
twice, independently, by first setting the initial value of 5a to zero whilst leaving 5(p to take its 
Bunch-Davies vacuum form: 


54>, 5a 
5(j)',5a' 



( 20 ) 

( 21 ) 


and then vice versa. We normalize the number of e-foldings to be iV = 0 at horizon exit, and so 
plot our results around this value. 

The code is split into four sections, each solved successively with the end values to each one 
used as the initial conditions in the next: 
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1. Inflation: This covers the period from N = 0 through to when the inflaton crosses the 
minimum of the potential, at which point we switch the decay, on. 

2. Inflaton decay: Covering the period through the first part of reheating, but before the 
secondary field has begun to decay. 

3. Overlap decay: The inflaton still contributes a significant amount to the overall energy 
density but the secondary field, a too has started to decay, so we switch on T„. 

4. Secondary field decay: Finally, we switch off the evolution of the (f> field altogether as it 
is so difficnlt and time consuming to follow the vastly different scales involved in both this 
and much smaller secondary field, cr. This then continues until the power spectra settles on 
a specific value and all the energy density of the universe is held within the radiation. 

We take the decay parameters to be T*^ = 10 “®Mpl and T®" = 10 “^^Mpl, and the masses of the 
scalar fields to be = 10 “^Mpl and ruo- = lO^^^Mpp. These values are chosen to be close to 
those in Ref. |26j . and compatible with the limits in Eq. (11) of Ref. [38] 

While the first two sections can be integrated in a matter of minutes, the latter sections can 
take some considerable time to track the oscillatory phases throughout decay of the secondary 
field. This is due to the small scales involved in comparison to the inflationary phases. Since the 
sudden decay approximation does not hold in this case |39|, the later sections are crucial in our 
numerical procedure in order to obtain an accurate result. Even in the standard case, with a = 0, 
we find that the ratio of a to the other components in the universe, rjeo evaluated when H = T^ 
is not the true point at which rdec reaches its maximum value. An improved (and increased) value 
can be attained slightly before this at Ff ~ 3 — ST*^ or, more accurately still, read off from its 
numerical maximum. The need for following the decay in full becomes even more apparent when 
we look at the results, in the next section. We find that oscillations in the non-adiabatic pressure 
perturbation towards the end of curvaton decay play an important role too. 

IV. RESULTS 

A. The case: CTmin = crini 

For this section we set cTmin = Cini hr order to exclude any evolution of the secondary field 
during inflation (see Figure [^. By keeping cTyain = CTini = 0.1 and comparing to a standard 
curvaton scenario, for which a = 0 and rdec — 0.18 (Figure [^, we see a significant change in the 
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Figure 1: Top: The evolution of the relative energy density in each species, for both the a = 0 
and a = —0.005 cases, which overlap throughout. Bottom left: The background evolution of a for 
both a values and for cTmin = = 0.1. Bottom right: The evolution of the effective Planck mass 

for a = —0.005 when U min = oiai = 0.1. 


amplitude of the final power spectrum for a given k. In Figure we plot the final twelve efolds 
as the inflaton decays, followed by the radiation dominated and secondary field dominated phases 
in terms of both (the dimensionless power spectrum) and hP^ad- In this and later plots in the 
paper, we take k = 0.05 Mpc“^. This shows the influence of the non-adiabatic pressure on the final 
power spectra; the SPaad survives for around an efold longer and has a maximum amplitude of up 
to roughly 100 times that of the standard case. Figure shows that this increase in amplitude 
is not due to a more dominant secondary field, as the value of flo- at the time of decay remains 
roughly constant (the change is of the order of 0.1%). 

For the case with Uini = 0.1, as above, we see in Figurej^that the dependence on a is independent 
of sign. This might be expected due to the boost in power spectrum coming as the secondary field 
oscillates and decays. Terms such as / and /, which both contribute to hP^ad, effectively average 
out as their sign changes back and forth. For this reason, in the case of cTmin = we shall now 
only look at the effect of |a|. For the cases when iTmin 7^ crinij which we will consider in the next 
section, this may no longer remain true as the sign of a will introduce a scaling of the final field 
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Figure 2: The power spectrum of the curvature perturbation (left) for both a = 0 and 
a = —0.005 cases and the associated (5Pnad (right). 



Figure 3: The amplitude of the power spectrum as a function of a normalised to the a = 0 power 

spectrum: Vc^(a)/Vc_(Q). 

value, (Tend) which can in turn affect the final power spectrum. We also observe a slight dip on 
either side of o; = 0, for which the amplitude decreases before increasing again. We have checked 
that this is not a numerical artefact. We do not have a physical explanation for this dip and the 
complexity of the governing equations makes it difficult to address this question analytically. 

Finally, in these simple cases, it is useful to compare how (Tend affects the outcome for various 
values of |a|. We will specifically focus on (Tend = 0.1,0.2,0.3 which results in rjec values of 
0.17,0.45 and 0.62 respectively. For these final values of a we would expect a varying increase in 
the amplitude of the power spectrum arising simply from the standard curvaton results. At the end 
of inflation we find = 3.01 x 10“^^ and this value is boosted by factors of 2.51,11.13 and 13.85 
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Figure 4: Amplitude of the power spectrum as a function of a for three different values of 

Uend = 0.1,0.2,0.3 

respectively by the end of curvaton decay with a = 0. In Figure]^ this small boost is apparent in 
the values at a = 0 but is insignificant in comparison to the subsequent amplitude increases as we 
increase a from 0. The results for each Cend diverge for increasing a due to the fact that for each 
(Tend we also have Umin = CTend, SO that the difference between the true minimum (cr ~ 0) and the 
local minimum associated with /(cimin) increases as the values for cjend increase. 


B. The case: cTmin ^ cjini 

In the more general case we have two possibilities, namely a > 0 and a < 0. This choice plays a 
role in the evolution of a during inflation which can in turn affect (Tend and rjec- By choosing a to 
be negative, we can pull a towards its local minimum before decay; a positive a has the opposite 
effect and pushes it away. This second case soon becomes unworkable for values of a approaching 
0.05 or greater due to the exponential increase apparent in the value of a. Due to the symmetry of 
OL shown in Figure]^ this need not be too concerning, however, as once a reaches its final value at 
the end of inflation we can still study the subsequent effects purely by using negative values. The 
only benefit of using a > 0 comes in the ability to further vary the trajectory of a to test that our 
results are independent of it. This can be done by setting cr to a number of different initial values 
and using a to control its final value in order to compare results. 

In each case we find that the power spectrum is dominated by the value of (Tmin with a lesser 
but still noticeable dependence on (Tend- This is most simply demonstrated by Figurein which we 
show two cases, both with a^ia = 0.1 but with ccini = 0.1 and 0.3 respectively. We let a run over 
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Figure 5: Left: The power spectra for varying a with a\n\ = 0.3 and fJmin = 0.1 (red) in 
comparison to the case of crini = (T min =0.1 (blue). Right The background trajectories for a for 

each of these cases. 


the same values used previously which, for Uini = 0.3, gives various values of (Tend- 0 > a > —0.03 
results in 0.3 > Uend > 0.1, while a < —0.03 gives Uend = 0.1 (see the right hand side of Figure]^. 
It is clear from the left hand side of Figure that while a remains small, the final amplitude of the 
power spectrum differs from that of the Uini = cTmm case. This can be explained by the observation 
that in these cases the field has not had enough time to reach cjend = 0.1 due to the smallness 
of a. For larger a, however, the two cases converge because Uend now equals 0.1 in each of these 
examples. 


C. The case: cTmin = 0 


Finally, we give an example which demonstrates both that the evolution of a throughout in¬ 
flation has little to no impact (other than the dependence on Uend) on the final power spectra 
amplitudes and that we get no noticeable boost when (Tmin = 0. Here we take iTini = 0.05 and 
a = {0.0056, 0.011,0.0145} which gives Uend = {0.1, 0.2, 0.3}, respectively. From these we hnd that 
the amplitude of the hnal power spectrum is multiplied by factors of 1.03,1.29 and 1.51. However, 
in comparison to the factors involved in the standard curvaton-like case given in Section IV A for 
these values of Uend we see that the changes represent additional increases of only 1 — 3%. These are 
insignificant when taking into account the usual approximations inherent in the curvaton model 
and those increases found earlier in the paper for Umin / 0. 
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V. CONCLUSION 

In this paper we have studied a model of the early universe consisting of two scalar fields: the 
inflaton and a second field which controls the stabilisation of the Planck mass. We work in the 
Jordan frame, for which the second field is non-minimally coupled to gravity; this choice allows us 
to deal with the decay of the fields in the usual way. We have investigated numerically the effects of 
this coupling on the power spectrum of primordial fluctuations. It has previously been shown that 
a non-minimally coupled scalar field can induce changes in the curvature perturbation on super¬ 
horizon scales via the introduction of terms proportional to /, / and / in the non-adiabatic pressure 
perturbation, JPnad [H]- Here, we have quantified this effect. We have shown that it can play an 
important role on the amplitude of the power spectrum in a non-minimally coupled curvaton-like 
case, in which the secondary field decays only after inflation is complete. Allowing the effective 
Planck mass to evolve in such a way, even by the smallest of amounts, leads to dramatic changes 
in the amplitude of the final power spectrum in comparison to the standard curvaton scenario. 

The effect of this amplitude boost can also be linked to the spectral index, and tensor-scalar 
ratio, tts by parameterising the power spectrum as m 

= (1 + R)v'f^ , (22) 

where 

-pfo-) 

' c 

This gives 

4e — 2r]^ , 16e , , 

Us - I = -2e + 2rj„ - 1 ^ = TTR ’ 

using the usual definitions of the slow roll parameters, evaluated at horizon crossing. The important 
point to note here is that n* and depend only on the ratio, R, not the mechanism by which 
the curvaton, or curvaton-like field, sources the final curvature perturbation. This places tight 
constraints on the values that a, cTmin and hence f{cr) can take according to the latest Planck 
data mi- We soon find ourselves with a spectral index approaching 0.98 - as in the pure curvaton 
limit - as R becomes large with relatively small changes in /. This is also largely independent of 
any evolution in a during inflation because the inflaton dominates the universe at horizon crossing, 
when both the slow roll parameters are evaluated and the tensor perturbations freeze in. 
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It remains to be seen whether a similar scenario to the one discussed will arise from fundamental 
theories of particle physics. If so, it will have an impact on inflationary model building in such 
theories. 
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